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Abstract 

We discuss the possible existence of an observational window, in the TeV- 
PeV energy range, for the detection of prompt neutrinos from the decay of 
charmed particles produced in cosmic ray interactions with the atmosphere. 
We calculate the event rates for muon and tau neutrinos of heavy quark, 
mostly charm, origin. We argue that their prompt fluxes are observable in 
a kilometer-scale neutrino telescope, even though the calculations are sub- 
jected to large uncertainties, which we quantify. We raise the possibility that 
a small component of prompt neutrinos may already be present in the ob- 
served samples of current experiments. We also discuss the interplay of the 
predicted fluxes with those produced by the flavor oscillation of conventional 
atmospheric neutrinos, and by anticipated cosmic sources. 
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The quest to search for sources of cosmic neutrinos beyond the sun, to search for the 
particles that constitute the cold dark matter, and to exploit other science opportunities 
ranging from astronomy to particle physics, led to the commissioning of large volume high- 
energy neutrino telescopes. While AMANDA and Baikal are taking data similar and 
much larger instruments are contemplated [§-0. These detectors collect the Cherenkov 
radiation emitted by charged secondaries (electrons, muons and taus) produced in neutrino 
charged current interactions in the ice or water surrounding the optical sensors. In order to 
filter the background of atmospheric muons created by cosmic-ray interactions, only upward- 
going neutrinos which traverse the Earth are monitored. The first mission of a new neutrino 
telescope is to calibrate the detector on the known flux of atmospheric neutrinos J/J. Up to 
about 10 TeV, the main source of atmospheric neutrinos is the decay of pions and kaons in 
the atmosphere produced in the interactions of cosmic rays with the Earth's atmosphere; we 
will refer to them as constituting the "conventional" atmospheric neutrino flux. At higher 
energies, these mesons interact rather than decaying into a neutrino because of the increasing 
lifetime of the parent mesons. Therefore the semileptonic decay of very-short lived charmed 
particles becomes the dominant atmospheric source, giving rise to the "prompt" neutrino 
flux. The energy dependence of prompt neutrinos follows the cosmic ray spectrum whereas 
the spectrum of conventional neutrinos is steeper by one power in energy because of the 
competition of decay with interaction of the parent particles. The prompt neutrino flux is 
independent of zenith angle whereas conventional neutrinos are preferentially produced in 
the rarified atmosphere at large zenith. 

In this letter we discuss the possible existence of an observational window for prompt 
neutrinos in the TeV-PeV energy range. We calculate neutrino induced muon and tau event 
rates, recognizing that the fluxes are subject to large uncertainties arising from the combi- 
nation of extreme atmospheric cascade parameters with different charm production models. 
We will argue that the prompt fluxes are observable in high-energy neutrino telescopes and 
analyse the interplay of the predicted fluxes with those produced by the flavor oscillation of 
conventional atmospheric neutrinos, and by anticipated cosmic sources. 

Prompt neutrinos are produced by production and semileptonic decay of heavy quarks 
produced in cascades initiated by cosmic rays in the atmosphere. In order to appreciate 
the main ingredients in calculating the production of prompt leptons in the atmosphere, 
we write an approximate solution to the cascade equations ||, valid for energies below the 
charm critical energy (f« 10 8 GeV): 

HE) = Z pi (E) Z a (E) W °)» C 1 ) 

where <&i(E) is the flux of secondary leptons (I = fi, v e and v^) from the decay of charmed 
hadrons (i = D,D S or A c ) produced by cosmic ray hadrons. The charm production 
spectrum- weighted moment Z pi (E) describes the inclusive charm production cross section; 
Zn(E) is the charm decay spectrum- weighted moment which represents the semileptonic 
three-body decay kinematics; A P (E) and \ P (E) are, respectively, the proton attenuation 
and interaction lengths in air ||. The primary cosmic ray flux at the top of the atmosphere 
& P (E,0) is assumed to be composed of protons although the validity of the superposition 
approximation guarantees that the calculation extends to heavier nuclei. Above the charm 
critical energy (the energy where interaction and decay of the charm particle compete in 
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the evolution of the cascade), the charm interaction length Xi(E) must also be taken into 
account and the complete solution has to be considered. What the prompt tau-neutrino 
flux (/ = u T ) is concerned, it suffices to calculate a generalized decay Z-moment Zn(E), in 
Eq. fll]) which accounts for the decay chain: D s — > ru T followed by r — > u T X (where X is 
either a fiv^ pair or a meson) fTOfl . 

Detailed investigation of the charm induced muon-neutrino flux j XI] reveals the large 
uncertainties in the prompt component resulting from the imprecise knowledge of atmo- 
spheric particle showering parameters: primary cosmic ray spectrum $ p , spectrum- weighted 
moments Z iu , attenuation and interaction lengths, A p , X p and Aj. Further uncertainty is 
associated with the extrapolation to high energy of a variety of models describing the ac- 
celerator data on charm production such as the Quark Gluon String Model (QGSM) |i"2| , 
Recombination Quark Parton Model (RQPM) |Ljj and perturbative QCD (pQCD) |RJ. In 
the end, we conclude that the calculated fluxes, for both and is T , are subjected to uncer- 
tainties of up to one order of magnitude for a given charm model. When comparison is made 
of alternative charm production models, the spread reaches two orders of magnitude. We 
will therefore define an allowed range between the maximum and minimum prompt neutrino 
fluxes. 

There are two contrasting approaches in calculating charm production cross sections. 
The first is to blindly apply perturbative QCD. This is a questionable procedure because, 
at the extreme energies considered here, large logarithms connected with the small ratio 
of the charm quark mass and the square root of the center of mass energy, x = m c /s/s, 
are likely to modify the calculation. For instance, this procedure, even though subject to 
large ambiguities, routinely underestimates the production of strange particles at accelerator 
energies. It may therefore be more appropriate to calculate charm production in the context 
of QCD-inspired models that have been patterned to accommodate accelerator data on 
strange particle production. QGSM and RQPM fall into this category and predict similar 
and larger extrapolations for the charm production cross section at higher energies. 

The detection rates of upward-going prompt neutrinos, i.e. the rates of neutrino induced 
muons and taus, per year per effective detection area (in 2ir sr) , can be estimated from our 
flux calculations following the prescription of Gandhi et al. |15| . We first concentrate on the 
operating neutrino telescopes, BAIKAL NT-200 g] and AMANDA-II 0, and subsequently 
consider the proposed km 3 experiment, ICECUBE |}J; see Table |. Calculations use the 
maximum and minimum allowed prompt fluxes. We show for comparison the expected rates 
for muon-neutrinos from conventional atmospheric decays. Absorption in the Earth of the 
highest energy neutrinos is taken into account. From Table p] we anticipate a situation 
where more then a thousand prompt muon-neutrinos above 1 TeV may be extracted from 
ICECUBE data using their flat zenith angle distribution as a signature. Furthermore, we 
predict about 200 prompt tau-neutrino candidates which also produce characteristic showers 
in the detector ||16|| . Above 10 TeV, the 20 remaining tau-neutrinos will have no counterpart 
from the tau-neutrinos resulting from the oscillation of conventional atmospheric neutrinos, 
a source of which we consider further on. 

Our calculations raise the possibility that prompt neutrinos may be observed by high- 
energy neutrino telescopes, especially because models extrapolating to the larger fluxes are, 
at least in our opinion, more robust. We next discuss the interplay of the predicted fluxes 
with those produced by the flavor oscillation of conventional atmospheric neutrinos, and by 
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anticipated sources of cosmic neutrinos. 

In Figure [l] we present prompt fluxes corresponding to the allowed range between 
maximum and minimum predictions |Tl[ , along with the conventional atmospheric flux p4 



Among the potential candidates of other high-energy neutrino sources [|17]], we will con- 
centrate on cosmic accelerators such as active galactic nuclei and gamma ray bursts which 
produce the highest energy photons observed, and are theorized to produce the highest 
energy cosmic rays. Their anticipated diffuse neutrino flux is bounded by experimental in- 
formation on measured high-energy fluxes of cosmic rays, specifically neutrons, and gamma 
rays. The two upper curves in Figure |TJ labeled MPR (solid and dotted lines), represent the 
allowed range of neutrino fluxes from sources that are optically thick or thin to the emis- 
sion of neutrons |18| , respectively. The two intermediate straight lines labeled WB (dashed 
and long-dashed), represent bounds imposed on sources which produce the highest energy 
cosmic rays [|l9j. The higher flux allows for cosmological evolution of the sources. 

We note that the energy where the prompt neutrino flux exceeds the conventional one 
is roughly 10 TeV. Separation of the two components can be made, taking advantage of the 
fact that the conventional atmospheric flux has a characteristic and calculable zenith angle 
distribution, while the prompt component is isotropic. 

If cosmic sources exist near the MPR bounds, the prompt muon-neutrino window is 
definitely closed. Neutrino telescopes will detect cosmic neutrinos at very high rates and the 
prompt component is a background to be dealt with. The AMANDA experiment is already 
exploring this range of fluxes 0. Between 10 5 GeV and 10 6 GeV the MPR bound for thin 
sources has essentially the same slope as the prompt flux. Since the bound is an upper limit, 
it may happen that the cosmic sources actually produce as many detectable neutrinos as 
atmospheric interactions do, and it will be difficult to disentangle both components. The 
situation is quite different when considering the extragalactic flux subjected to the WB- 
bounds. For evolving sources, there is room for an excess of prompt neutrinos in the range 
20-300 TeV. For non-evolving sources scenario, then the prompt window may be open up to 
about 2 PeV. 

The allowed band for atmospheric charm-induced tau-neutrinos is presented in Figure |2|. 
In the absence of neutrino flavor oscillations, atmospheric z/ r 's from charm are the only 
source of tau-neutrinos and the prompt window is completely open. We remind the reader 
that cosmic beam dumps produce negligible fluxes of tau-neutrinos. Recent measurements 
of the Super-Kamiokande experiment |2(J do however favor non-vanishing neutrino masses 
and produce compelling evidence for — > v T flavor transitions, with maximum mixing angle 
(sin 2 29 = 1) and a square mass difference of Am 2 = 3.2 x 10~ 3 eV 2 . We therefore consider a 
scenario where the conventional atmospheric flux in Figure [T] is modified by this oscillation 
mechanism. In the calculation we consider the range log 10 (Am 2 /eV 2 ) = 2.5 ±0.5. Similarly, 
the bounds bracketing the extragalactic diffuse muon-neutrino flux can be translated into 
estimates of bounds on a cosmic tau-neutrino flux by adopting a typical value for the ratio 



of v T to vacuum flavor conversion, over astronomical baselines, of F Ut / u = 0.5 pi 
The resulting fluxes are also shown in Figure The appearance of a tau-neutrino below 
1 TeV is a clear indication of atmospheric neutrino oscillations from a conventional 
flux. Above ~ 2 TeV, we face again the interplay between extragalactic bounds and the 
atmospheric prompt component. For MPR scenarios, the prompt window may be closed; 
for WB-bounds the prompt tau-neutrino appearance may dominate extragalactic oscillation 
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fluxes up to 50 TeV (for evolving sources), or even 500 TeV (for non-evolving sources). 

In conclusion, we observe that the calculations of prompt neutrino fluxes are subjected 
to large uncertainties, arising from possible extreme cascade parameters combined with 
different charm production models. We identified the possibility that prompt fluxes are 
observable in a kilometer-scale neutrino telescope, particularly in the region from 2 TeV 
to 2 PeV. At the upper range of the predictions, a small component of prompt neutrinos 
may already be present in the observed sample of high-energy atmospheric neutrinos of the 
AMANDA experiment 0. 

Detection of the atmospheric prompt-^ will provide unique information on heavy quark 
interactions at energies not accessible to particle accelerators. The better characterization 
of the charm induced neutrino component will increase the discrimination power of neutrino 
detectors at higher energies. Knowledge of atmospheric charm production is essential in 
identifying the origin of a detected v T above a few TeV, whose appearance may also result 
from neutrino flavor oscillations. 
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TABLES 



TABLE I. Upward-going muon and tau event rates per year arising from charged current 
interactions of atmospheric neutrinos in ice or water, for different neutrino telescopes' effective 
areas and thresholds. Prompt flux calculation based on maximum (MAX) and minimum (MIN) 
range limits described in the text. 
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BAIKAL NT-200 


1 TeV 


2 x 10 3 m 2 


22 


3 











AMANDA-II 


1 TeV 


3 x 10 4 m 2 


330 


44 


2 


7 





ICECUBE 


1 TeV 


1 km 2 


11000 


1470 


53 


216 


11 


ICECUBE 


10 TeV 


1 km 2 


170 


157 


3 


22 


1 
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FIGURES 




FIG. 1. Comparison of several contributions to the high-energy z/^ flux. The band represents 
the allowed range between maximum and minimum atmospheric prompt neutrino fluxes. Thick 
dashed line is the conventional atmospheric flux. Solid, dotted, dashed and long-dashed lines 
correspond to upper bounds imposed to diffuse extragalactic neutrino flux by the observation of 
high-energy cosmic-ray and gamma-ray spectra, according to different scenarios explained in the 
text. Fluxes are multiplied by E 2 . 
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FIG. 2. Comparison of several possible contributions to the high-energy v T flux. The band rep- 
resents the allowed range between maximum and minimum atmospheric prompt neutrino fluxes. 
Thick dashed lines are the result of maximum mixing flavor oscillation of the conventional atmo- 
spheric Ufj, flux, for Am 2 around the Super-Kamiokande value (thick dot-dashed line), as explained 
in the text. Solid, dotted, dashed and long-dashed lines correspond to the upper bounds in Fig- 
ure 1, subjected to vacuum flavor oscillation, averaged in transit to Earth. Fluxes are multiplied 
by E 2 . 
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